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We study field theoretical models for cosmic strings with flat directions in curved space-time. 
More precisely, we consider minimal models with semilocal, axionic and tachyonic strings, respec- 
tively. In flat space-time, isolated static and straight cosmic strings solutions of these models have 
a flat direction, i.e., a uniparametric family of configurations with the same energy exists which is 
associated with a zero mode. We prove that this zero mode survives coupling to gravity, and study 
the role of the flat direction when coupling the string to gravity. Even though the total energy of the 
solution is the same, and thus the global properties of the family of solutions remains unchanged, the 
energy density, and therefore the gravitational properties, are different. The local structure of the 
solutions depends strongly on the value of the parameter describing the flat direction; for example, 
for a supermassive string, the value of the free parameter can determine the size of the space-time. 
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^ i T. INTRODUCTION 

F ^ . In several field theoretical models which allow for string-like objects it has been observed that families of solutions 
which have different local behaviour of the fields appear, but possess the same value of the energy per unit length. 
These zero modes do not disappear when coupling the models minimally to gravity. Therefore, despite the fact that 
the global properties of the family of solutions will remain unchanged, the gravitational properties might change, 
since, after all, gravity is intrinsically locally defined. 

The simplest model that is frequently used to describe cosmic strings is the U(l) Abelian-Higgs model. The 
| models we study in this paper are somewhat more complicated than this one, but will be closely related to it since 
the solutions can (sometimes) be seen as embedded Abelian-Higgs strings. Semilocal strings are solutions of an 
SU{2) g iobai x U(l)i OC ai model which - in fact - corresponds to the Standard Model of Particle physics in the limit 
sin 2 9 W = l^where 9 W is the Weinberg angle. The simplest semilocal string solution is an embedded Abelian-Higgs 
!— 1 ' solution [J Q ■ A detailed analysis of the stability of these embedded solutions has shown [H, H[ that they are unstable 
(stable) if the Higgs boson mass is larger (smaller) than the gauge boson mass. In the case of equality of the two 
masses, the solutions fulfill a Bogomolnyi-Prasad-Sommerfield (BPS) [j| bound such that their energy per unit length 
', is directly proportional to the winding number. Interestingly, it has been observed [!, H| that in this BPS limit, it 
is possible to find a one-parameter family of single static and straight cosmic string solutions: the Goldstone field 
I/"") | can form a non-vanishing condensate inside the string core and the energy per unit length is independent of this 
• value, which itself is related to the width of the string. These solutions are also sometimes denominated "skyrmions" 
and have been related to the zero mode associated with the width of the semilocal strings present in the BPS limit. 
Semilocal strings have been studied in cosmological settings both in the context of their formation [H-IH, network 
CNI ' properties [l(| and their CMB implications [Tl| . 

The second model that we are investigating in this paper is a field theoretical model for so-called axionic and 
tachyonic strings. This model, introduced in Blanco-Pillado et al. [l2j], originally describes an unstable D-brane anti- 
D-brane pair and allows for the existence of non-singular BPS strings. In fact, depending on the boundary conditions 
imposed on the matter fields, the model accommodates three different types of strings 13]: 0-strings (tachyonic), 
s-strings (axionic), and a hybrid of both types. As in the semilocal model, each of these types of strings has a family 
of solutions of the same energy parametrized by one parameter. For tachyonic and axionic strings, respectively, this 
parameter is associated with the "width" of the strings - very similar to what happens in the semilocal model. On the 
other hand, for hybrid strings it measures the contribution of the tachyonic string in relation to the axionic string. 
These strings, in particular the tachyonic strings, share some qualitative features with semilocal strings. 
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The gravitational properties of field theoretical cosmic string solutions have also been discussed in the literature. 
The most thorough study has been done for Abelian-Higgs strings by minimally coupling the Abelian-Higgs model to 
gravity [13, EE)]. Far away from the core of the string, the space-time has a deficit angle, i.e. corresponds to Minkowski 
space-time minus a wedge [l6j . The deficit angle is proportional to the energy per unit length of the string. If the 
vacuum expectation value (vev) of the Higgs field is sufficiently large (corresponding to very heavy strings that have 
formed at energies bigger than the GUT scale), the deficit angle becomes larger than 2tt. These solutions are the so- 
called "supermassive strings" studied in [I?], IH[ and possess a singularity at a maximal value of the radial coordinate, 
at which the angular part of the metric vanishes. Gravitating semilocal strings have first been studied in fl9j . 

In this paper we reinvestigate the gravitating semilocal strings and point out further details. Moreover, we in- 
vestigate gravitating axionic and tachyonic strings focusing in all cases on the role of the zero mode. Our paper is 
organised as follows: in Section [Til we describe the models and give the equations of motion. In order to simplify 
the calculations we will use a standard result in the context of supergravity theories: cosmic string configurations 
preserving a fraction of the supersymmetries in a supergravity theory saturate a BPS bound, and the corresponding 
Einstein equations admit a first integral, the so called gravitino equation, (see for example |20j). The models we will 
discuss here can be embedded in J\f — I supergravity, [l2L l2ll ] , and moreover the corresponding cosmic string solutions 
leave unbroken half of the supersymmetries of the original theory. Here, following the results in [22j, we will rederive 
the first integrals of the Einstein equations for the present models treating them as non-supersymmetric theories, and 
therefore, without making any explicit reference to supersymmetry. The rest of the field equations will be obtained 
using a BPS-type of argument valid for static cylindrically symmetric configurations in a gravitational theory, that is, 
searching for the conditions required to minimize an appropriate energy functional. The resulting equations will turn 
out to be a set of first order differential equations for both models analog to those found in flat space-time @, [l2j , 
and in particular we recover the results already found in for the semilocal case. 

In section IIIII we discuss the numerical results and conclude in Section IIVI 



II. THE MODELS AND THEIR FLAT DIRECTIONS 



The models we are studying have an action of the following form 

S = --Mp J d i x yS-dct gR + J d^x -Act g C m , (1) 

where the gravitational strength coupling is given in term of the reduced Plank mass M~ 2 = 8irG. The space-time 
metric has signature (— , +, +, +). The Ricci tensor scalar follows conventions from [20j . 

We will study two types of matter Lagrangians £ m , the semilocal model [l], 0, and the axionic D-term model 
presented in First we will discuss the gravity part of the model and then we will describe the matter content 
in the subsections below. In particular, we will point out what type of string-like solutions these models can 
accommodate and discuss the role that gravity plays in connection to the zero modes. 

In the present work we will focus on static cylindrically symmetric configurations invariant under boosts along the 
axis of symmetry which, without loss of generality, we can take to be the z— axis. The most general line element 
consistent with these symmetries is: 

ds 2 = -N 2 (r)dt 2 + dr 2 + L 2 (r)d<p 2 + N 2 (r)dz 2 , (2) 

where (r, ip, z) are cylindrical coordinates. The non- vanishing components of the Ricci tensor R 1 ^ then read 

ff o _ (LNN'Y 2N" L" (N 2 LJ 

where the prime denotes the derivative with respect to r. With our conventions the Einstein equations have the 
following form 

R,», = -M- 2 (T ia ,-±g la ,T\ , (4) 
where T = T£ is the trace of the energy-momentum tensor which is given by 

T »» = ~ 2 £)gT^ + ■ ( 5 ) 
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In [22[ it was shown that the energy momentum tensor of the semilocal model has a very simple form for field 
configurations satisfying the BPS bound. In particular, the only non-vanishing components of T£ are proportional to 
the field strength of an auxiliary vector field , which is a function of the gauge boson and the scalar fields involved 
in the string configuration: 

T; = T; = 0, Tl = T z z =±L-\d r A*-d v AZ), (6) 
As we shall see below this statement is also true for the axionic D-term string model proposed in [l^ . 

If the energy-momentum tensor is of the form ([6]) the metric function N(r) becomes a constant [ID, l24j - and thus 
choosing the coordinates conveniently, we can always set it to one, N{r) = 1. To obtain the equation of motion for 
the remaining metric function L(r) we use the (^-component of (J3J) which reads 

^ = -M- 2 (t* - \t} = hd- 2 T . (7) 

In order for the metric to be regular at the z— axis we need to impose the following boundary conditions 

L(0) = , L'(0) = 1 . (8) 

For cylindrically symmetric field configurations it is possible to choose a gauge where the vector field A^ satisfies 
A® = and A]? = A^ (r) and therefore, as long as the BPS bound is saturated, we can find a first integral to the 
Einstein equations 1 : 

L" = ±M~ 2 (A*y =► L' = 1±M- 2 A%. (9) 

Here the integration constant is fixed by the boundary conditions (0, since regularity also requires that ^4^(0) = 0. 
This first integral of the Einstein equations allows us to express the deficit angle of the string configurations in terms 
of the asymptotic value of the vector field A^ [23[ 

6 = 2tt(1 - L'\ r=00 ) = t2ttM- 2 A%\ r=ao . (10) 

In order to discuss the existence of zero-modes in these theories we need an appropriate definition of the energy which 
is valid for general curved space-times. We use the same definition as in [20J, which is valid for time independent 
configurations, and is obtained adding a Gibbons-Hawking term to the action: 

E = -S - S GH , S gh = M 2 V [ y/-detg K. (11) 

JdM 

Here K is the trace of the second fundamental form of the metric at the space-time boundary dM, and g a p is 
the metric induced at dM. Note that field configurations minimizing this energy functional also extremize the action, 
and therefore are solutions to the equations of motion [25| . In order to obtain a finite value for the energy we calculate 
the integrals in (TlTj) only over the plane orthogonal to the string, i.e. over the r and ip coordinates. Then, choosing 
the boundary to be a cylinder centered on the z— axis with an arbitrary large radius r — ¥ oo we have (after setting 
N(r) = 1) 



s/-detgK = L'(r), (12) 
and thus the Gibbons-Hawking term can be written in terms of the deficit angle 

Sgh - Ml f ^-detg K = 2n (L'| r=00 - L'\ r=0 ) - -A/ p 2 S. (13) 

JdM 

Due to the Einstein equations, the action S vanishes on shell; and we find the following relation between the energy 
of the string and its deficit angle: 

E = Mld = T2nA^\ r=00 , (14) 



1 In the context of supergravity theories this equation is know as the gravitino equation [20 
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where the second equality can be obtained using equation (JTUJ). In the following subsections we will arrive to the 
same result following a Bogomolnyi type of argument Q which involves rewriting the energy functional as a sum of 
positive terms plus a boundary term. 

It is possible to show that for field configurations with the same symmetries as the ones we discuss here the definition 
of the energy (fTTj) agrees with the one used in [19| . 



E = J d 2 x y/-det gT$ . (15) 

In order to check this is sufficient to note that using the ansatz for the metric ([2]), the Einstein-Hilbert term becomes a 
boundary term which is exactly cancelled by Sqh, and that any contribution to the action involving time derivatives 
must be zero. 



A. Semilocal strings 



First, we will reconsider the semilocal model which possesses a S , J7(2) g i b a i x £^(l)iocal symmetry 0, Q and is given 
by the following matter Lagrangian density 



(16) 



where F^ v = d^A^ — d^A^ is the field strength tensor of the U(l) gauge field and £> M $ = (e?p — iA^)<& is the 
covariant derivative of the complex scalar field doublet $ = (0i, </> 2 ) T . The constant g denotes the gauge coupling, A 
the self-coupling of the scalar fields and \/£ is their vacuum expectation value. 

If the couplings satisfy the Bogomolnyi limit g = and working in the gauge Aq — A z = 0, the energy functional 
(1111) for the stationary string with a translational symmetry along the z— axis can be given by 



E = J drd^L(r) (D r ^ D r <S> + L~ 2 D V <P + ^L.- 2 F rtp F rv + y (£ - * f $)' 

which can be rearranged in the following way: 



(17) 



E = 



drdipL(r) (D r $ ± LL- 1 D( P $)T(D r $ ± iL^D^) 



T / drckpFZ 



(18) 



where we have used the explicit form of the metric ((5]) after setting N(r) = 1, and we have introduced the auxiliary 
vector field A^ and its field strength 2 : 



A* = U&D^-Dp&ty-tA,. 



F^^d^-d^A*. 



(19) 



In order to obtain finite energy configurations (1171) we have to require that far away from the center of the string the 
fields are in the vacuum, and that the covariant derivatives of the fields vanish 



$ + $| r ^oc=£, D^$| r _,oo=0. (20) 
Therefore, using Stokes' theorem, we find the following lower bound for the energy of the cosmic string 



E>T drd(p F° = ±£ / drdip F, 



(21) 



where the last integral is the total magnetic flux trapped inside the string. Any field configuration which saturates 
this bound is at a local minimum of the energy functional and thus, as we argued in the previous subsection, it must 



2 In the context of A/ - = 1 supergravity, is known as the gravitino U(l) connection. 
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be a solution to the equations of motion. In particular, the bound is saturated by cosmic string configurations which 
satisfy the BPS equations: 

D r $ ± iL^D^ = 0, L- 1 F r9 T9 2 = 0. (22) 
We use the following ansatz for a static straight cosmic string lying along the z— axis [H, 0]: 

0i = \TU{ry ntp , 02 = V?Mr)e™^ , A v = ±«(r), (23) 

with all the other components of the gauge field set to zero. Here n and m are the winding numbers of the fields <f>\ 
and 4>2 respectively and, without loss of generality, we will assume 3 \n\ > \m\. By using the rescalings r — > r/(g-\/^) 
and L(r) — > L(r) / '(ff\/£) the BPS equations read 

/' + ^^/ = 0, h' + ^^h = 0, ^ + (/ 2 + /, 2 -l) = 0, (24) 

which have to be solved subject to appropriate boundary conditions that result from imposing (1201) and regularity at 
the origin 0, Q 

fL + hl = l, u(r -> oo) = |n| , (25) 
/(0) = , v(0) = , C=o(0) = or /i m#0 (0) = , (26) 

where /oo = f(r = oo) and ft-oo = ft.(r = oo). The choice of signs of the winding numbers n and to in (|24p ensures 
that these conditions can be met. Indeed, the BPS equations (|2"2")) with the upper sign can only be solved provided 
n,m > 0, while the lower sign requires n, to < 0. 

The first two BPS equations (12"4")) imply that the profile functions / and h must be related to each other [l], [|J 

logfc = log/-(|n|-M)/;^ + * h = c-f-exp((\m\-\n\)J^j (27) 

and therefore, there is a one-parameter family of solutions characterized by a real constant n = log c . As in flat space, 
this family is degenerate in energy, i.e. the total energy does not depend on the parameter k. This can be checked 
inserting the ansatz for the gauge boson (|2"3"]) and its boundary conditions in (|21[l : 

E = ±£ [dipA 9 \ r=00 = 27r|n|£ . (28) 

For field configurations satisfying the BPS equations (f2"2"j) it is easy to check that the energy momentum tensor takes 
the form (0, with given by (fl^l) . Then, the profile function of the metric L(r) can be obtained from ^j, which 
using the ansatz ((23]) takes the form 

L' = 1 - a 2 ((|n| - «)/ 2 + (\m\ - v)h 2 + v) , (29) 

where a = M~ V? i s the vacuum expectation value of ^ measured in Planck masses. From the previous two equations 
it is possible to recover the relation (fl"4"]) between the deficit angle of the BPS strings and their energy 

6 = M- 2 E = 2ir\n\a 2 . (30) 

Note that the choice of signs made above for the winding numbers ensures that the energy, and thus the deficit angle, 
are positive. 



B. Axionic and tachyonic strings 

In this subsection we are going to consider the cosmic string solutions of the axionic D—term model studied in 
fiat space-time in |12l Il3l |. The model describes the dynamics of two complex scalar fields, the tachyon </>, and an 
axio-dilaton S = s + ia (s > 0), coupled to a U(l) gauge field A^. The lagrangian density reads 



Ci = -D^D»<t> ~ RggD^SD^S - -^F^F^ - 9 - (£ + 28K S - q^f (31) 



3 The case n = m can always be rotated into a Nielsen-Olesen string using a SU(2) transformation, and therefore will not discuss it. 
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where F^ v — F M „ = d^A v — d v A^ is the U(l) field strength, and the covariant derivatives of the tachyon and the 
axio-dilaton are given by D^(f) = (<9 P — iqA^cj) and D^S = d^S + i2SA^ respectively. The model depends on three 
continuous parameters: the gauge coupling g, the charge of the axio-dilaton 5, and £ which determines the expectation 
value of the fields. The constant q is an integer which represents the U(l) charge of the tachyon. The lagrangian 
density also involves the derivatives of the Kahler potential associated with the axio-dilaton field K(S, S), which is a 
real function chosen to have the following form: 

K (S, S) = -M 2 \og{S + S). (32) 

The quantities K$ and K S g denote the derivatives of the Kahler potential with respect to the fields S and S, and 
they are given by 

K s = -Ml— , K sq = Ml- . (33) 

p S + S ' ss p (S + S) 2 v 1 

This lagrangian density is the bosonic part of a supersymmetric model with a D— term scalar potential, and therefore 
the couplings satisfy the Bogomolnyi limit [26| . Thus, without imposing further constraints in the couplings, the 
energy of a static configuration 

E = J drd^L{r) (^D r 4>D r ^ + L^D^D^ + K sS D^SD»S + ^L,- 2 F rv F rv + 9 — (f + 25K S - q^f\ , (34) 

can be written in the Bogomolnyi form as follows 

E = Jdrd(pL(r) \D r <j> ± iL^D^ 2 + K sS \D r S ± iL~ l D v S\ 2 

+ ^(L- 1 F rvT g 2 (Z + 2SK s -q\cj } \ 2 )) 2 } T J drdipF^ , (35) 

where we have defined the composite vector field and its field strength 

iM 2 

AZ^itfD^-tDpfi-^siDfJS-D^-ZAp, F*=d u Af-d„A2. (36) 

As in the semilocal model, if we restrict ourselves to finite energy configurations, we have to impose that far away 
from the string core the fields are in the vacuum, and that the covariant derivatives vanish 

( g |0| 2 - 2MCs) r _> 00 = £ , £ M 0Uoo , D^Slr^co = . (37) 

With these boundary conditions we find a lower bound for the energy similar to (|2f p . and the corresponding BPS 
equations (which ensure that the bound is saturated) read 

{D r ±iL- 1 D v )<f> = 0, (D r ±iL- 1 D v )S = 0, L^ 1 F rip =p g 2 {£ + 25K S - q\4>\ 2 ) = 0. (38) 

As we mentioned earlier in the paper, the energy momentum tensor reduces to a very simple form for field configura- 
tions which satisfy the previous equations. Indeed, the only non-vanishing components are given by 

T l t = T* = ±L- X F^ (39) 

and therefore as we discussed before the Einstein equations admit the first integral ©• 

In order to solve the BPS equations we use the ansatz for the matter fields proposed in [l2j], which represents a 
static straight cosmic string along the z— axis: 

4>=VUqf(r)e in * s~\r) = £/(<5M 2 ) h{rf 

a = 25rrup A v — ±v(r) . (40) 

After rescaling the radial coordinate and the metric profile function, r — > ?"/(<?vD> L — > L/(g^), the BPS equations 
become: 

/' + ^H/ = 0, ft'+o^^O, V - + (f + h*-l) = 0, (41) 
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where a — M p 1 \Ji/q- The signs of the winding numbers n and m are fixed by requiring /(r) to be regular and 
h(r) 2 > for r — > 0. As in the semilocal case, we can use the first two equations to find a relation between the 
tachyon and the dilaton field 

1 C rlr 

= 2(\n\-q\m\) / -— -2tog/ + «, (42) 
(ah) J L(r) 

and therefore the solutions are parametrized by an arbitrary constant k. Actually, this model admits three different 
families of cosmic string solutions depending on the boundary conditions at r —¥ oo: 

• ^-strings (tachyonic) 

In this type of strings the magnetic flux trapped in the core is induced by the winding of the tachyon field, which 
must satisfy \n\ > q\m\ in order to solve the BPS equations. The profile functions have the following asymptotic 
behavior 

f(r -> oo) -> 1 , h(r -> oo) -> , v(r -> oo) -> \n\/q (43) 

so that the tachyon field acquires a non- vanishing expectation value far away from the core while the function 
h(r) tends to zero. For the solutions to be regular we also have to impose the following boundary conditions at 
the core of the strings r — > 0: 

/(0)=0, w(0)=0, /4 =0 (0)=0, or ft m#o (0)=0. (44) 

In the present work we will only consider the case m = 0, where the profile function of the axion-dilaton, h(r), 
can be non-zero at the core creating a "condensate" . Actually, as in the flat space-time analysis [l2| , the value 
of h(r = 0) will be a free parameter related to «, which determines the width of the string. This family of 
solutions is degenerate in energy 

E = ±i J dipA v \ r=00 = 27r|n|~ , (45) 

and therefore the zero-mode survives the coupling to gravity. For simplicity, in our numerical analysis we will 
set the constant 5=1. 

• s-strings (axionic) 

In this case the magnetic flux inside the strings is induced by the winding of the axio-dilaton, S. The behaviour 
at infinity is 

f{r — > oo) — > , h(r — !• oo) — !• 1 , v(r — >• oo) — >• \m\ . (46) 

It is now the dilaton which acquires a non-zero expectation value far from the core, while the tachyonic field 
tends to zero far from the core. Thus, the role of the tachyonic and the dilatonic field is exchanged. These 
strings are solutions to the BPS equations provided |n| < q \m\. Regularity at the origin imposes the following 
boundary conditions 

h(0)=0, w(0)=0, /; =0 (0)=0, or / n#0 (0) = 0. (47) 

In the next section we discuss the properties of this family of solutions for the case n = 0, where the value of 
the profile function of the tachyon at the center of the string, f(r = 0), is a free parameter related to K. Again, 
this family of solutions is degenerate in energy 

E = ± tJ dipA v \ r=00 = 27r|m|£, (48) 

and thus the zero-mode still exists after coupling the model to gravity. Note that s— strings have a tension q 
times larger than 0— strings. As for the tachyonic strings, we will restrict the analysis to the case q = 1. 

• Hybrid-strings 

In this case, both the tachyon and the dilaton field contribute to cancel the scalar potential far away from the 
core, i.e., they both acquire a finite vacuum expectation value for r —¥ oo. This can only happen provided the 
windings satisfy \n\ — q \m\. Thus, in this case we have the following boundary conditions at infinity 

fl + hl = l, »( r ->oo) = M = M, (49) 
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where f(r —> oo) = foo and h(r — > oo) = hoo. Provided the previous constraints are satisfied, the value of /oo 
is a free parameter, which will be again related to k. In this case k is not a measure of the width of the string, 
instead, it is related to the relative contributions of the tachyonic and the axionic string to the tension of the 
string . The corresponding boundary conditions at the core of the string are given by 

/(0) = 0, «(0)=0, h(0)=0. (50) 

The zeromode associated with the parameter k is not normalizable and thus we will not discuss it any further in 
the present work. Indeed, if we promote the parameter k to be time dependent the corresponding effective action 
for the zeromode gets a quadratically divergent contribution A 2 , where A is a cutoff which, in a cosmological 
setting, could be given by the distance to the closest cosmic string. 

Finally, the form of the equation for the metric profile function ^ can be found using the ansatz (HU1) : 

L' = 1 - a 2 [(|n| - qv)f + q(\m\ - v)h 2 + qv] . (51) 

Similarly to the semilocal case, from the boundary conditions discussed above, it is straightforward to find a relation 
between the deficit angle of these strings and their energy: 6 = M~ 2 q E. 



III. NUMERICAL RESULTS 

We have solved the BPS equations for both models, subject to appropriate boundary conditions, focusing on the 
interplay between the flat direction and gravity, in order words, we studied how the local field configuration changed 
the energy density and the properties of the solutions, even though the global energy was unchanged. 

A. Semilocal strings 

The gravitating properties of semilocal strings were first studied in 19]. Besides completing that analysis, we 
include this case in the present work because it shares several properties with the other cases to be studied. 




FIG. 1: Profiles of a typical n = 1 semilocal string solution uncoupled to gravity, with different values of the free parameter: 
h(0) = 0.1 (left) and h(0) = 0.9 (right). 

Fig. [T] shows two typical n = 1 cosmic string solutions to this model, when the coupling to gravity is not included. 
The metric function L = r in this case. These figures illustrate the effect of varying the value of the condensate h 
at the center of the string r = 0, i.e. the free parameter of the family of degenerate solutions. For larger values of 
h{r = 0) the profile functions reach their asymptotic values farther away from the center of the string and thus, as 
we anticipated in the previous subsection, the width of the string core increases. 

The form of the profiles when coupling to gravity, can be seen in Fig. [5J Close to the core, the profile functions 
have the same form as in a Minkowski background: 

I(r)«rt..., /M«/ rl n '+... , h(r) « ho r |m| + . . . , v(r) « \r 2 + . . . . (52) 



FIG. 2: The profiles of a gravitating n = 1 semilocal string (a = 0.5) with the value of the free parameter h(0) =0.1 (left) and 
h(0) = 0.9 (right). 



However, far away from the core, r — > oo, the metric takes the form 

ds 2 « -dt 2 + dr 2 + (1 - \n\a 2 ) 2 r 2 dip 2 + dz 2 , 



(53) 



which corresponds to a conical space-time with a deficit angle 5 — 2Tr\n\a 2 , and the behavior of the profile functions 
in the same limit is given by the following asymptotic expansion: 



fir) 



1- 



1 



3(|m|-M) 

l-l„.l^ 



h(r) 



v(r) ss |n| — (|n| — |m|) 



> (54) 



where ro is a parameter related to ho which determines the width of the string core. Setting the constant a = we 
recover the corresponding expansions in flat space. 






h(0)=0.1 




h(0)=0.5 




h(0)=0.9 







FIG. 3: The function L for different configurations. On the left, h(0) = 0.5 while we vary the values of a; on the right a = 0.5 
while we vary the value of h(0). 

The figures show that the metric function L does indeed depend on the value of the condensate. We can see 
the dependence more clearly in the right plot of Fig. [3] where we have displayed the metric function L for different 
values of the condensate keeping constant the value of a. In the left plot of Fig. [3] we show the profile function L 
in different cases where we keep the value of the condensate constant and vary the value of the a parameter. Note 
that the deficit angle (|10p increases with increasing a, and also that as a approaches zero, the space-time becomes 
Minkowski, where L(r) = r. 

One way of understanding this phenomenon is by looking at the energy density for each one of these configurations, 
as shown in Fig. |H where we only vary the value of the condensate. Even though the total energy (the area below 
these lines) is the same, the energy density is different, and affects the metric function with different strength. Indeed, 
since the width of the string core increases with the value of the condensate, the energy density spreads and the metric 
profile function L reaches its asymptotic behavior farther away from the string. Note also that the deficit angle for 
each of the configurations is exactly the same, and therefore the slope of the profile function L in the limit r — > oo 
should be the same regardless of the value of the condensate. 

For BPS cosmic strings the energy density is closely related to the magnetic field F rip 



(55) 
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FIG. 4: Energy density for different values of the free parameter h(0), with fized a = 0.5. 



and thus from ([5"2"jl and (|5l|) it is trivial to find the limiting form of the energy density 



£{r -> 0) w | 



£(r — > oo) 



2g(H-H) 2 
(i 



n cr 



2(|m|-|"l) 

ro 



(56) 



A dramatic effect happens when considering supermassive strings. In this case the string is massive enough to make 
the function L turn round and become zero at some finite value of r* (see Fig [3] for the case with a — 1.02), in other 
words, the deficit angle S becomes larger than 2ir. The value of the condensate in the core determines the extent 
to which this solution exists; it decides the "size" of the universe. Once again, the total energy does not change, 
i.e., the integral for the energy density curve from r = up to the other point r* where L{r*) — is independent 
of the value of h(r = 0). In Fig. [5] we have depicted the values of r* for a fixed coupling constant with respect to 
the value of the condensate at the core h(r = 0). As the value of the condensate increases the string core (where 
space-time is approximately Minkowski) becomes wider. In the limiting case where h(r = 0) — > 1 the space-time 
becomes Minkowski everywhere and the point r* tends to infinity. 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

h(0) 



FIG. 5: Values r* at which the function L = for a given coupling constant (a = 1.02) with respect to the value of the 
condensate at the core h(0). 



B. Tachyonic strings or cj> strings 

These solutions share many properties with the semilocal strings, the main difference being the factor of h 3 in the 
equation (|4"Tj) instead of the h in equation This translates into the function h tending to zero logarithmically. 

As in the semilocal case, the tachyonic field / is responsible for the formation of the string, whereas h is responsible 
for the condensate. Close to the core of the string the profile functions behave as in flat space-time [l2j|, in particular 

/(r)«/ rl"l+... , h(r)- 2 ^ h 2 -2a 2 q\m\logr+... , v(r) « \r 2 + . . . . (57) 

The value h(0) is a free parameter, related to the string width, which leaves the total energy unchanged and modifies 
the field configurations slightly. A typical n — 1 cosmic string configuration can be seen in Fig. [SJ which shows that, 
as in the case of semilocal strings, for larger values of the condensate the core width increases. Note also how the h 
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FIG. 6: Profiles of a gravitating tachyonic string (a = 0.5) with the value of the free parameter h(0) =0.1 (left) and h(0) = 0.9 
(right). 



field tends to zero slowly. Actually, all the profile functions approach their values at infinity logarithmically, as can 
be seen from their asymptotic expansion 



(l-|n|a 2 ) 



2a 2 (|n| — q|m|) logr 



+ ..., /(r)«l- 



(l-|n|a 2 ) 



4a 2 (|n| — q|m|) logr 



. , v(r) fa 



(1- \n\a 2 f 



q Aa 2 q(\n\ — q\m\) log r 



(58) 

The term appea ring in the numerator of the three expressions (1 — |n|a 2 ) was not obtained in the flat space-time 
analysis done in [121 ] , since this is a consequence of the asymptotic form of the space-time metric 



From 



ds 2 ps -dt 2 + dr 2 + (1 - \n\a 2 ) 2 r 2 d(p 2 + dz 2 . 
we can immediately extract the form of the energy density far away from the core: 



£{r) 



M 2 (l-\n\a 2 ) 
2(|n| — g|m|)r 2 log 3 r 



(59) 



(60) 



h(0).0.1 - 
h(0).0.5 
h(0).0.9 - 



FIG. 7: Energy density for different values of the free parameter h(0) with fixed a = 0.5. 

Most of our results for the semilocal strings can also be obtained here. The global energy, and thus the deficit 
angle, are independent of the value of the h field at r = 0, implying that the slope of L remains unchanged for large 
values of r; however, as the energy density spreads due to the presence of the condensate, the behaviour of the metric 
profile function L changes close to the core, and in particular reaches its asymptotic behavior farther away from the 
string center (see Fig. [5]). 

Once again, the deficit angle increases for larger values of a, and in the case of supermassive strings, the metric field 
develops a zero far from the core, making the space-time closed. Fig. O depicts the points r* where the L develops a 
second zero for different values of the condensate at the core, keeping the value of a fixed. 



C. Axionic strings or s-strings 



For type of strings, the roles played by the tachyon field / and the dilaton field h are interchanged, the latter being 
responsible for the formation of the strings, and the former giving a measure of the width of the string (see Fig. I10p . 
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FIG. 8: The function L for different models: left, h(0) = 0.5 and various a; right, a = 0.5 and various h(0). 




FIG. 9: Values of the points r* at which the function L — for a given coupling constant a — 1.1 with respect to the value of 
the condensate at the core h(0). 



Close to the core, the approximate form of the profile functions is also given by (|57l) . but in this case q\m\ > |n|, and 
in particular Fig. HU1 corresponds to the case n — 0. Here /(0) plays the role of the free parameter which fixes the 
width of the string. The following asymptotic expansion shows the behavior of the profile functions in the opposite 
limit, r — > oo 



2(|r.|-g|m|) , , 2(|n|-g|m|) 



r \ l-^lmla 12 \ f T \ l-^lmla 13 | Tl | — q | TO | / T 



\roJ 2 \r J qa z \r 0/ 

_ (61) 

This result is slightly different than the one obtained in a Minkowski background [12], as can be seen from the 
appearance of the factor (1 — g|m|a) in the exponents. As for the case of the tachyonic strings, this correction is 
related to the form of the space-time metric far away from the core of the string: 

ds 2 w -dt 2 + dr 2 + (1 - q\m\a 2 ) 2 r 2 dip 2 + dz 2 . (62) 

From (|61[) we obtain the asymptotic expansion of the energy density far away from the core: 

T2 



2M 2 {\n\ - q\m\) 2 f r 



£ W = hi I 2 \ 2 2 - + • • ■ • ( 63 ) 

(1 — q\m\a z ) z r z \ro J 

The energy density of these strings is close to that of the semilocal or tachyonic strings, in that there is a concen- 
tration of energy next to the core, although at much shorter distances than in those cases. This effect has its origin 
in the factor Re(S') -2 of the kinetic terms of the axio-dilaton field ([Mf. which is responsible for the divergence of 
Re(S) at the center of the s— strings. As in the two types of strings previously discussed, the details of the shape of 
the energy density depend on the value of the condensate. 

Once again, the global energy of these configurations does not change with respect to the value of the condensate 
at the core, the zero mode survives. In this case, solutions for supermassive strings that have a closed space-time 
associated with them can also be obtained, and the dependency of the values r* at which the closing occurs can be 
found in Fig. QTJ 
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r r 

FIG. 10: The profiles of a gravitating axionic string (a = 0.5) with the value of the free parameter /(0) = 0.1 (left) and 
/(0) = 0.9 (right). 




o 1 1 1 1 1 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

HP) 

FIG. 11: Values of the points r* at which the function L = for a given coupling constant a = 1.1 with respect to the value 
of the condensate at the core /(0). 

IV. CONCLUSIONS AND DISCUSSION 

In this paper we have studied field theoretical models for cosmic strings with flat directions in curved space-time. 
Specifically we have considered the effects of the gravity coupling on models with semilocal, axionic and tachyonic 
strings, which were known to have flat directions in flat space-time. In this work we have focused on solutions of a 
single static cylindrically symmetric string with invariance under boosts along the axis of symmetry. Although the 
models studied are very different in character, some of the solutions obtained are very similar. 

The cosmic string solutions have been found using a BPS-type of argument, which consists in minimizing an 
energy functional appropriate for static cylindrically symmetric field configurations in a gravitational theory. The 
field configurations saturate the corresponding BPS bound provided the fields satisfy a set of first order differential 
equations, similar to those found for Minkowski space. Following [22j we have shown that, if these conditions are met, 
then the energy momentum tensor acquires a very simple form, which allows to find a first integral to the Einstein 
equations. This result is also interesting in the context of supergravity theories, where this first integral is known as 
the gravitino equation, and it is known to exist whenever the cosmic string configuration preserves a fraction of the 
supersymmetries. However, since the derivation we have presented here makes no reference to supersymmetry, our 
results show that cosmic string solutions saturating a BPS bound can satisfy a gravitino-type of equation, regardless 
of the fraction of supersymmetries of the model that is broken by the string. 

As in the case of the non-gravitational version of these three types of strings, the set of differential equations 
characterizing them admit a one-parameter family of solutions which are degenerate in energy. Therefore, we have 
proved that the zero mode survives the coupling to gravity in the three cases. We have found numerical solutions to 
the BPS equations for the three types of strings, paying special attention to the space-time metric, and in particular 
we have characterized its dependence on the value of the family-parameter and the string tension. 

Even though the free parameter does not change the total energy, and therefore, neither the deficit angle of the 
string, it changes the shape of the energy density. This is due to the fact that the free parameter changes the 
profiles of the fields forming the string, in particular it changes the width of the string core, which in turns modifies 
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the energy distribution. As gravity depends on the energy density, not only on the global value of it, the energy 
distribution changes the metric function L which, the wider the string, the farther from the string center it reaches 
its asymptotic behavior. Thus, we show that different values of the zero mode do change the metric properties, even 
though they do not change the global characteristics, such as the deficit angle. This effect becomes very apparent 
when the string is massive enough to have a deficit angle larger than 2tt. In those cases the L function turns round 
and becomes zero at some finite value of the radial coordinate, rendering the spatial directions transverse to the 
string closed. As the free parameter associated with the zero mode changes the shape of the L function, its value 
acts as a modulus which fixes the size of the compact space-time dimension. 

We would like to mention possible generalizations of this work. On the one hand, it would be interesting to study 
whether the zero-mode associated with the string width would survive the coupling to gravity in general situations, i.e. 
situations involving several strings with a random motion. It is known that the semilocal model admits static parallel 
muti- vortex solutions with flat directions which survive the coupling to gravity [23 |. and the low energy dynamics of 
the corresponding zero-modes has been studied in flat space-time [27|- Given the similarities between the semilocal 
strings and the axionic and tachyonic strings we expect that performing similar analyses for the later cases is possible, 
and in particular to find static multi-vortex solutions with flat directions even after including the gravity coupling. 
On the other hand, the study of the models in this paper was performed statically. The analysis of the dynamics of 
these models would also be very interesting, for example to investigate whether the zero modes could be excited, or 
whether the zero modes could actually be dynamical. 
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